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We investigate strongly interacting dense matter and neutron stars using a flavor-SU(3) 
approach based on a non-linear realization of chiral symmetry as well as a hadronic fiavor- 
SU(2) parity-doublet model. We study chiral symmetry restoration and the equation 
of state of stellar matter and determine neutron star properties using different sets of 
degrees of freedom. Finally, we include quarks in the model approach. We show the 
resulting phase diagram as well as hybrid star solutions for this model. 
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1. Introduction 

The study of strongly interacting matter at extreme conditions of temperature and 
density is at the forefront of modern nuclear physics. This regime covers the physics 
of ultrarelativistic heavy-ion collisions as well as important aspects of nuclear as- 
trophysics. Whereas the determination of the phase structure of excited strongly 
interacting matter at high temperature and baryon densities, as pursued in heavy- 
ion collisions, involves relatively high temperatures, the complementary study of 
the structure of compact stars is directly related to the properties of very dense and 
rather cold matter. 

In order to study all these regimes in a unified model approach we developed 
an effective chiral SU(3) model that can be studied over the whole relevant range 
of chemical potentials and temperatures. Reflecting the correct degrees of freedom 
of the system after the transition to a quark-gluon plasma state we include quarks 
and the Polyakov loop as order parameter for the deconfinement phase transition. 
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2 S. Schramm and V. Dexheimer 

2. The Hadronic Model 

The model used in om calculations is an extended a—u chiral model, based on a non- 
linear realization of chiral symmetry, that includes the lowest SU(3) multiplets of 
baryons and mesons as well as an effective field that mimics the QCD scale anomaly. 
The expectation values of the scalar isoscalar fields correspond to the non-strange 
isoscalar (cr field), isovector (6) and strange (ss) chiral quark condensates (()• 

The interactions between baryons and the scalar and vector mesons are linear 
and read (assuming static systems, including only the time component of the vector 
fields) 



where w, p, and (f) are the non-strange isoscalar, isovector and strange vector fields, 
respectively. The index i sums over the baryon octet (N, A, E, S). Additional terms 
include mass terms and quartic self-interactions of the vector mesons, whereas the 
self-interactions of the scalar mesons induce the spontaneous breaking of chiral 
symmetry. The effect of non-zero current quark masses is modelled by introducing 
an explicit chiral-symmetry breaking term (for more details on the full Lagrangian 
see . The vacuum masses of the baryons are generated through their coupling to 
the chiral condensates. 

For the baryon-vector couplings giuj and pure /-type coupling is assumed 
as discussed inHJ, gi^ = (n^ — n^q)gs , gi<t> = "('^s ~ "Dv^Ss' > where g^ denotes 
the vector coupling of the baryon octet and the number of constituent quarks 
of species i in a given hadron. The resulting relative couplings are in accordance 
with additive quark model constraints, coupling only the vector mesons to strange 
baryons. 

All parameters of the model are fixed by either symmetry relations, hadronic 
vacuum observables or nuclear matter saturation properties (see^. In addition, 
the model also provides a satisfactory d escr iption of realistic (finite-size and isospin 
asymmetric) nuclei and neutron stars Using this approach the equations of 
motion are solved for isospin symmetric nuclear matter as well as star matter by in- 
cluding leptons and requiring charge neutrality. Fig. 1 shows the resulting equation 
of state of cold symmetric nuclear and star matter. The parameters were slightly 
tuned from the fit values used in calculations of finite nuclei to explore the range 
of possible maximum star masses, while retaining a reasonable description of sat- 
urated nuclear matter. The nuclear compressibility has a still acceptable value of 
297 MeV The results for three choices of baryonic degrees of freedom are shown, 
which comprise (a) nucleons, (b) nucleons and hyperons, and (c) nucleons, hyperons 
and spin 3/2 baryonic resonances. As expected, with increasing number of particle 
species the equation of state becomes softer at larger densities. Using these equa- 
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tions of state in a calculation of the properties of a spherical and static neutron 
star, integrating the Tolman-Oppenheimer-Volkoff (TOV) equations we obtain 
maximum star masses between 1.93 and 2.12 solar masses depending on the num- 
ber of degrees of freedom taken into account (Fig. 2). One can observe that even 
including A baryons in the calculation rather large star masses with typical radii of 
about 12km are possible. Looking into the corresponding particle densities, shown 
in Fig. 3, the hyperon content of the star is relatively small with A baryons ap- 
pearing at about 3 times nuclear matter groundstate density and the at 4 po, 
respectively. The latter result is due to the fact that in the model the S potential 
in nuclear matter is positive (Us — 5.35MeV). In case one includes the baryonic 
spin 3/2 decuplet, that is A baryons for the densities occurring in compact stars, 
the particle mix looks very different as can be seen in Fig. 4. The A~ is populated 
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Fig. 1. Equation of State of isospin-symmetric nuclear matter and stellar matter including nu- 
cleons, nucleons and hyperons, and in addition spin3/2 baryons. 
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Fig. 2. Mass-Radius diagram for neutron star solutions using different EOS as shown in Fig. 1. 
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above 2 po and the hyperons are shifted to larger densities. Note, however, that the 
results depend strongly on the largely unknown value of the vector coupling of the 
A to the Lo meson, here taken to be identical to the nucleon-oj coupling. 

One possible effect that might occur in neutron stars is meson, in particular 
K~ , condensation ^. Using our model, in a calculation of the critical densities 
for the onset of the condensation in neutron and proto-neutron stars we find that 
for cold stars the condensate sets in at values of about 6 times nuclear matter 
saturation density whereas for proto-neutron stars the value is even higher. These 
values are beyond the maximum central densities inside of the stars such that kaon 
condensation does not affect stellar properties I^EI. 
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Fig. 3. Particle abundances as function of density for a star including nucleons and hyperons. 
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Fig. 4. Particle abundances as function of density for a star including nucleons, hyperons, and A 
resonances. 
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3. The Parity-Doublet Model 

Using an alternative approach of formulating a chirally symmetric Lagrangian we 
consider the parity doublet model originally discussed in Here, one adopts the 
so-called "mirror assignment" for the positive and negative parity nucleon states 
(iV_l_ and N_, we restrict ourselves to SU(2), in this case), grouping both states in 
the same multiplet. Under left- and right-handed transformations SUl{2) x SU{2)r 
transformations the two nucleon fields ■01 and "02 transform non-trivially as: 

i>lR > R-lplR , V'li * LlpiL , (3) 

'4'2R > Ltp2R , 'lp2L >■ R'ip2L ■ (4) 

In a similar way as a and tt are parity partners in the linear realization of chiral 
symmetry this allows for a chirally invariant nucleonic mass term in the Lagrangian: 

'71o(V'275V'l - V'l75V'2) = rno{lp2L'iplR - 1p2R.'iplL ~ 1plL'ip2R + '4>1R^iI^2l) , (5) 

where toq represents a bare mass parameter. After diagonalizing the quadratic terms 
in the Lagrangian one obtains the unmixed fields, and -/V_, that are the standard 
nucleon and its parity partner. In the limit of chiral symmetry restoration (cr = 0) 
both nucleonic states are degenerate, but attain a finite mass mo. This doubling of 
degenerate nucleonic states is the characteristic feature of this type of chiral model. 

Using the parity-doublet ansatz combined with a linear sigma model including 
the w meson (for achieving the correct nuclear matter saturation properties) and the 
p meson (for reproducing the phenomenological value for the asymmetry energy of 
32 MeV) one can determine the corresponding field values and the equation of s tate 
by solving the equations of motion of the model in mean-field approximation 12121 
In Fig. 5 the resulting a field, the scalar condensate, is shown. The field becomes 
smaller with increasing density tending towards chiral symmetry restoration. The 
structure in the curve for star matter at about 2 times po corresponds to the onset 
of populating the parity partner of the neutron. This can directly been seen in Fig. 
6. Above 4 po the parity partner of the proton shows up. 

Inserting the equation of state into the star calculation one obtains a maximum 
star mass of 1.85 solar masses. Including sizeable vector meson self-interaction terms 
leads to very low star masses in contradiction to values of observed stars. Going be- 
yond the mean-field approximation we repeated the calculations using the Relativis- 
tic Hartree Approximation (RHA), which in general generates comparable results 
slightly lowering the nuclear matter compressibility '■2- 

In general, for the parity partners to be populated in the interior of the star 
a relatively low vacuum mass of the N~ baryon below 1380 MeV is required. In 
addition, in order to maintain reasonably small values for the compressibility one 
has to assume a high value of the mass parameter niQ > SOOMeV. A full analysis 
of the question whether with these values one can still be in accordance with low- 
energy hadronic data like pion- nucleon scattering is still an open question 
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4. Including Quarks 

Lattice QCD calculations show that the phase transition from hadrons to quarks 
and gluons is a smooth crossover at vanishing chemical potential. Some lattice 
calculations at finite chemical potential /is, which are still notoriously difficult to 
perform, sugges t a critical end-point of a line of first-order phase transitions in 
the T-fj,B plane Calculations that connect hadrons and quarks in two different 
model approaches are not able to reproduce such a transition behavior. Therefore 
we extend our hadronic m odel by i ncluding quarks and an effective field for 
the Polyakov loop, following I ^ I ^ I 3 1^ in order to describe the deconfinement phase 
transition in a single unified approach ^1^. The potential for the Polyakov loop reads: 

U = (aoT* + aifi% + azrVD^"^ + 03^0* In (1 - + 8$^ - 3$^). (6) 

The parameters are obtained by fitting the values to lattice QCD results at 
vanishing chemical potential. The /XB -dep endent terms are fixed by reproducing 
the critical end-point as computed in H^l xhe switch between quark and baryon 
degrees of freedom is achieved in an effective wav, shifting the baryon/quark masses 
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Fig. 5. Scalar condensate as function of density in the parity-doublet model. The mass for the 
parity partner N— is set to 1200 MeV. 
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Fig. 6. Densities of particle species as function of density in the parity-doublet model. 
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to high/small values for small/high Polyakov loop values using a simple additional 
mass term for quarks and baryons, such that their effective masses read: 

ml = gtaO- + gbSTsd + gtcC + mot + (7) 

ITT'q = ffg^Cr + gqST36 + gq^C + + (1 - , (8) 

where mg are small explicit mass terms and gh<i> and gq^ are the coupling constants 
of the Polyakov loop to baryons and quarks, respectively Using this approach 
we obtain the phase diagram shown in Fig. 7. Note that we also obtain a realistic 
first-order liquid-gas phase transition. The transition to quarks at zero temperature 
occurs at 4/9o- 

It is straightforward to calculate star masses in this model. The result in Fig. 
8 shows a maximum star mass of 1.92 solar masses. Using a Maxwell construction, 
assuming local charge neutrality of the star matter, the family of stable stars ends 
with the appearance of quarks and a significant softening of the equation of state. 
Assuming global charge neutrality one obtains a mixed phase of quarks and hadrons 
in the inner 2km core of the maximum mass star. 
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Fig. 7. Phase diagram in iig and T for the Hadron-Quark model. First-order transition lines for 
the chiral and deconfinement transitions as well as for the liquid-gas phase transitions are shown 
as calculated in the model. 
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Fig. 8. Mass-radius relation of compact stars in the Hadron-Quark model. The inset shows the 
different results for a Maxwell and Gibbs construction of the phase transition. The Gibbs con- 
struction generates a mixed hadron-quark phase that smooths the sharp transition. 



